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High-accuracy multicolor light curves of the binary system HD 189733, which contains an 
exoplanet, are analyzed. We have determined the radii of the star and the planet in the 
binary system, as well as the orbital inclination. The limb-darkening coefficients of the 
stellar disk were obtained in 10 filters in the wavelength interval AA = 5500 — 10500A The 
uncertainties of the fitted parameters were estimated using the differential-correction and 
confidence-area methods. The wavelength dependence of the limb-darkening coefficients is 
compared to the corresponding theoretical function for a model thin stellar atmosphere. We 
confirm the wavelength dependence of the exoplanet's radius found by Pont et al. [6j (at the 
la level). The exoplanet radius increases with decreasing wavelength, which seems to argue 
for the presence of an atmosphere around the planet. 



INTRODUCTION 

The results of fitting high- accuracy multicolor light curves for the star HD 209458, 
which is eclipsed by an exoplanet, were provided in pQ. There was a significant deviation 
between the observed and theoretical limb-darkening coefficients, even though the most 
conservative values were used to estimate the "external" uncertainties in the model pa- 
rameters [21 [3] . This confirms the discrepancy between observations and the theory of thin 
stellar atmospheres found in [I]. It is important to analyze limb darkening across other 
stellar disks that are eclipsed by their exoplanets to further study this effect. Pont et al. 
[SI E] obtained high-accuracy multicolor light curves of HD 189733, which is eclipsed by 
an exoplanet, and examined the spotted disk of the star. Pont et al. [5j [6] indicated that 
the exoplanet could have an atmosphere. The purpose of the current paper is to study the 
limb darkening of the star in the eclipsing system HD 189733 as a function of wavelength 
in detail using the data obtained in [5] . We analyze in detail the effect of uncertainties in 
the parameters determined using different methods to fit the eclipsing light curves, as in 

CP- 



FITTING METHOD 



The method used to fit transit light curves observed in eclipsing binary systems con- 
taining exoplanets is described in detail in pQ, and we will only consider the main concepts 
of the method here. We used a model consisting of two spherical stars in circular orbits, 
without any reflection effect or ellipsoidal effect. There are no spots on the surfaces of 
the stars. If the masses of a star and planet are m s = O.825M and m p = 1.15Mj up [7], 
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the mean relative radius of the Roche lobe for the planet calculated using the formula of 
Eggleton [16] is Rr/ci = 0.0531, where a is the radius of the relative orbit of the system. 

Since the relative radius of the planet is r p = R p /a ~ 0.0175 (see below), the planet fills 
its Roche lobe only to the degree /i p ~ 0.33, which is considerably less than 0.5. Therefore, 
a spherical approximation is quite satisfactory for the planet (neglecting flattening of the 
planet due to its axial rotation). The same can be said for the optical star. 

The amplitude of brightness changes due to the reflection effect in the optical should 
be less than 10 -5 magnitude at eclipse phases, which is negligible [1] . A weak wavelength 
dependance of the radius of the exoplanet in HD 189733 was suggested in [6] , which could 
be indicative that the planet has an atmosphere. The refraction of the stellar light in 
the exoplanets atmosphere can distort the light curve when the star is eclipsed by the 
exoplanet. The effects of light refraction in eclipsing binary systems were calculated by 
Kudzei [8]. These effects should result in small humps (brightenings) in the light curve 
before and after the eclipse, as well as in the middle of the eclipse. Since there are no 
such humps in the light curves of HD 189733, refraction effects in these light curves can 
be assumed to be negligible. 

When the star is eclipsed by its exoplanet, the light curve could be affected by grav- 
itational microlensing. Such effects were studied in [9]. Microlensing was shown to be 
significant only for the stars with exoplanets whose orbital sizes exceed 10 AU (with 
corresponding orbital periods of P > 10 yr). Therefore, our model with two spherical 
components is quite applicable to the light curves of HD 189733. 

We calculated the light curve using linear and quadratic limb- darkening laws to de- 
scribe the brightness distribution across the stellar disk: 



Here, p is the polar distance from the center of the stellar disk, Iq the brightness at 
the center of the disk, r s the radius of the stellar disk in radii of the relative orbit, and 
x the limb-darkening coefficient. The brightness at the center of the planetary disk, and 
accordingly the brightness at any point of this disk, are assumed to be zero. The planet 
eclipses the star when the orbital phase is 6 = tt. The orbit is taken to be circular, and 
the radius of the relative orbit a to be unity. The model does not include third light. The 
radius of the planet in relative orbital radii is denoted r p . The desired parameters of our 
model are the radius of the star r s , the radius of the planet r p , the inclination of the orbit 
to the plane of the sky i, and the limb- darkening coefficient x; if the limb darkening is 
described by the quadratic function, the parameters also include the coefficient y. 

The total light of the system is taken to be unity, and the observed brightness values to 
be normally distrubuted. The rms deviations of the observed brightnesses a are assumed 
to be known. 

To search for the optimal (central) values of the model parameters and generate their 
confidence intervals, we used a residual equal to the sum of the squared differences between 
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the observed and theoretical brightnesses (as functions of the desired parameters), divided 
by the dispersions of the observed brightnesses (the squares of the rms deviations a). If 
amodel fits the observational data, this residual is distributed according to a Xm law 
for the exact parameter values, where M is the number of observational points. This 
residual (which is a convex function of the parameters) achieves its minimum for certain 
parameters values r c s , rjj, i c , x c , y c , referred to as the central (or optimal) parameters. When 
the model fits the observational data, the distribution of the central parameter values can 
be considered to be normal in some vicinity of their mean values, if we neglect non- 
linearity of the parameter-dependent brightness in this neighbourhood. This is justified 
when a is small and the number of observational points M is large. 

The central parameter values are random values distributed normally in this ap- 
proximation. They can be used as a statistic in the differential-correction method or 
Monte-Carlo method (which, in a sense, is equivalent to the differential-correctionmethod 
PHI [EE])- The error intervals for the parameters are chosen to be centered on the optimal 
values of the corresponding parameters; their half-widths are equal to the mean square 
estimates of the rms deviations of the current parameters a es t from their central values, 
multiplied by a coefficient corresponding to the chosen confidence level 7 (la, 2a, 3a, etc.) 
[TT] . In this case, the model is considered to be perfectly valid; if we have more than 
one parameter, each individual parameter is characterized by its individual confidence 
interval, which encompasses the true value of that parameter with a given probability 7 
(independent of whether or not the true values of the other fitted parameters fall into their 
confidence intervals). We also assume that we can neglect non-linearity of the function de- 
scribing the dependence of the brightness on the parameters within approximately one rms 
deviation around the central parameters. Numerical experiments carried out for eclipsing 
systems show that the probability for all the true parameter values to simultaneously 
fall into their confidence intervals is a factor of 1.2 1.5 lower than the given confidence 
level 7 [10]. Since in the differential-correction method (or Monte-Carlo method) the 
obtained central parameter values are distributed statistically according to a normal law 
with the strict a priori assumption that the model is perfectly valid, and information on 
the dispersion of the observed brightnesses is not used here (the confidence intervals are 
constructed using the mean square estimates of the central parameters a 2 est rather than 
the dispersions of their central values), the uncertainties of the parameters in the differ- 
entialcorrection method (or Monte-Carlo method) are internal uncertainties. As a rule, 
these are substantially underestimated. Popper [J_2] indicated that this underestimation 
can reach a factor of three to five for eclipsing systems. 

In view of the above, we also used the confidencearea method to estimate the parameter 
uncertainties [13] . The advantage of this method is that the obtained confidence intervals 
for the parameters ensure that the exact values of all the parameters simultaneous fall 
into the corresponding confidence area (in a multidimensional parametric space) with the 
given probability 7. The probability that the exact value of an individual parameter falls 
into the corresponding confidence interval is greater than the given probability 7. In 
addition, the probability for the exact solution (a combination of all fitted parameters) 
to be covered simultaneously by all the confidence intervals is also greater than 7, since 
the volume of the parallelepiped in which the confidence area D is inscribed exceeds the 
volume of this area. 
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Therefore, the confidence-area method yields more conservative external estimates of 
the parameter uncertainties, which are a factor of a few greater than the error estimates 
obtained using the differentialcorrection method orMonte- Carlo method. 

We will use either a Xp or a Xm distribution to find the confidence area D JTUJ [TT] . 
We will use the difference between the residual distributed according to the x\a 1 &w (M 
is the number of observational points) and the minimal value achieved with the central 
parameter values as the Xp statistic, where P is the number of fitted parameters. If we 
neglect the non-linearity in the dependence of the brightness on the parameters (which is 
justified when there are many observational points M), these differences can be considered 
to be distributed according to the \ 2 p 1 &W f° r the exact parameter values, where P = 4 
with linear limb darkening and P = 5 with quadratic limb darkening. 

When using statistics distributed according to x\t or Xp> the confidence area D is taken 
to be a multidimensional set in the parameter space, for which the specified statistics are 
less than the quantile corresponding to the given confidence level 7. The probability for 
D to encompass the combined exact parameter values is then 7. The confidence intervals 
(uncertainties) for the parameters corresponding to a given confidence level 7 are the 
projections of the sides of the P-dimensional parallelepiped (in which the confidence area 
D is inscribed) onto the coordinate axes in the parameter space pU]. As was noted above, 
this method for constructing the confidence intervals ensures that they will simultaneously 
encompass the exact values of all the parameters with the given probability 7 (or higher) . 
Therefore, when specifying the confidence intervals (uncertainties) for the parameters, we 
should bear in mind that the adopted 7 (for the differential-correction or Monte- Carlo 
method) is equal to the probability for the exact value of each fixed parameter to be 
encompassed by the corresponding confidence interval; at the same time, the probability 
for the exact solution to be encompassed simultaneously in the all confidence intervals 
is less than 7. When we use the confidence area method, the probability for the exact 
value of a fixed parameter to be included in each individual confidence interval is greater 
than the adopted 7, whereas the probability that the confidence area D encompasses the 
combined exact values of all the parameters is 7. Therefore, the adopted confidence level 7 
specifies each of the individual confidence levels in the differential-correction method and 
Monte-Carlo method, and is related to the total confidence area D in the confidence-area 
method. 

Similar to the differential-correction method, using the x 2 p statistic assumes that a 
model is perfectly valid, and that the obtained asymptotic confidence area D never de- 
generates into an empty set. The confidence area D does not depend on the sample size, 
and the probability 7 can be specified for an area D that includes the exact solution (this 
probability is referred to as the confidence level). Since it is strictly assumed that the 
model is perfectly valid (the residuals are minimized over the parameters), the param- 
eter uncertainties obtained using the xp statistic can also be considered to be internal, 
though these are more conservative than the uncertainties obtained using the differential- 
correction or Monte-Carlo method. We must bear in mind, however, that the dispersions 
of the observational data are used to construct the x 2 p statistic; this is a significant differ- 
ence between the use of the xp distribution and the differential-correction method. This 
difference can also be qualitatively significant if the rms dispersion with unit weighting 
(distributed according to a reduced x 2 l aw ) differs significantly from unity; that is, when 
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there is a basis to think that the model does not fit the observational data well, or a 
strong correlation is suspected in the observational data. The parameter uncertainties 
obtained using the Xp distribution can also be considered external, since the dispersions 
of the external observational data are used here. 

The most conservative "external" parameter uncertainties are obtained using the Xm 
statistic, which does not assume the model to be perfectly valid and uses the dispersions 
of the observational data. It is also important that the same residual (the x\i statistic) is 
used both to search for the central parameter values and to find their confidence intervals. 
However, the confidence set D can be empty for certain significance levels a when using 
the Xm statistic. As a result of fitting the observations, only a less strong statement can 
be made concerning the true parameter values, since we did not assume by default that 
the model is perfectly valid when using the x\i statistic; namely, we can only indicate 
the probability (the significance level a) that wemistakenly reject the hypothesis that the 
combined true parameter values belong to a certain area D (rather than, strictly speaking, 
the probability 7 with which the confidence area D encompasses this combination of 
parameters). If we suppose that the model fits the observational data, the probability 
for the area D to encompass the exact values of all the parameters is 7 = 1 — a, where 
the probabilty 7 = 1 is referred to as the confidence level for the confidence area D. If, 
however, the model does not fit the observational data, this may mean that we make a 
second-order error in accepting the model (i.e., the model is false, but is accepted according 
to the statistical criterion). In this case, the equality 7 = 1 — a is only approximate for 
the confidence area. 

However, we must remember that, in most cases, we are dealing with a specific real- 
ization (a sample) of a light curve rather than the entire observed light curve (a general 
population). Therefore, difficulties that arise when trying to demonstrate the adequacy 
of a model may be associated with considerable statistical deviations of a specific obser- 
vational data, rather than shortcomings of the model itself. A light curve obtained at 
another epoch may match the model adequately, with no difficulties arising when fitting 
the observations. Therefore, strictly speaking, a specific realization used in fitting does 
not necessarily imply that the model is inadequate; it implies only that the available ob- 
servational data are not sufficient to ascertain whether the model is adequate or not. It is 
this that enables us to estimate the model parameters and their uncertainties even when 
the reduced x 2 1 X 2 e di * s significantly greater than unity. 

It follows from these considerations that searches for model parameters and their 
uncertainties in such problems should be accompanied by testing of the model adequacy. 
The Xm statistic can be used to test the adequacy of a model and verify how it fits the 
observational data (since the model is not assumed to be perfectly valid). Here, we must 
bear in mind that we accept the model according to a statistical criterion not because it 
is perfectly valid, but because there is no reason to reject it. 

We can judge how well the model fits the obesrvational data (whether the model is 
adequate enough) by estimating the significance level a = a starting from which the 
model can be rejected (the value is called the critical value of the significance level). 
The higher the critical value ao, the greater the probability that we make a first-order 
error when rejecting a model (that is, the weaker the basis for us to reject the model). The 
critical significance level «o is associated with the reduced x 2 f° r the minimal residual, 
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xled = XM M_p' n j which decreases with increasing a [TT] . 

Strictly speaking, the residual minimized over nonlinear parameters in nonlinear para- 
metric problems is not distributed according to a xii-p 1 &W ; an d only asymptotically ap- 
proaches this law as M — >• oo [13] . Therefore, the criterion for an adequate model that the 
Xred re d is close to unity can be used only when the number of observational points M is 
sufficiently large. It was indicated in [11] that, if xled = 1> the critical significance level is 
a ~ 50%, which corresponds to a very good model. In fact, we have the possibility in this 
case (when choosing a corresponding significance level a) to make up to 50% first-order 
errors in fitting if we reject the model; that is, we can err in every second case in rejecting 
the model. Hence, we have no reason to reject the model, and it can be accepted. When 
xl ed exceeds unity, the corresponding critical significance level is cto < 50%. In this case, 
we do not make a large number of firstorder errors in rejecting the model; i.e., we are 
correct in most cases when we reject the model, and we therefore have a basis to reject the 
model. If, however, we obtain Xred < 1; the corresponding significance level is «o > 50%. 
When the observational data are normally distributed, this situation is very unlikely, even 
for a model that is perfectly valid (see, for example, [IT]); therefore, the value Xred < 1j 
and the corresponding value a Q > 50%, can indicate the presence of a correlation in the 
input observational data, which may include not only random but also systematic errors. 

Finally, let us illustrate again why it is necessary to verify the model adequacy when 
deriving parameters from fits. An independent way of verifying the model adequacy 
follows from the fact that the number of observational points M usually greatly exceeds 
the number of parameters P, rule, the problem is highly overdetermined. Strictly 

speaking, one should have an exact light curve corresponding to the true parameter values 
to estimate their uncertainties. However, an exact light curve (corresponding to the 
perfect parameter values) is not available in reality; we have only an observed light curve, 
subject to observational uncertainties. The question thus arises of how well the observed 
light curvematches the perfect light curve. In other words, before finding the solution 
from a fit, we must answer the question of whether we have a sufficient basis to substitute 
the observed light curve for an exact light curve that is unknown a priori. 

We can answer this fundamentally important question after checking the models ad- 
equacy. If xled — 1' that is, a ~ 50%, there is no reason to reject the model, and the 
model can be accepted. An imagined perfectly exact light curve would optimally pass 
through the observational points of the real light curve, and we would have a firm basis 
to substitute the observed light curve for the perfectly exact light curve. The obtained 
values of the model parameters are close to their true values. Moreover, since the observed 
light curve optimally matches the perfectly exact curve and the contribution of systematic 
errors to the "exact" curve is negligible, we can obtain reliable uncertainty estimates that 
encompass the true parameter values with the given probability 7. We can estimate the 
parameter uncertainties using various methods, depending on the strictness with which 
we wish to judge the solution. These include the differential-correction or Monte-Carlo 
methods, which enable us to obtain "internal" uncertainties, the confidence-area method 
based on the Xm statistic, which yields "external" uncertainties, or the confidence- area 
method based on the Xp statistic (the intermediate case). 

However, if we obtain xled > 1 an d> correspondingly, a < 50%, we cannot make an 
unambiguous judgement about whether or not to accept the model. The modelmay prove 
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to be "bad" since we can adopt for the model only a significance level a for which the 
probability to make a first-order error is small; i.e., we must only rarely be mistaken 
when rejecting the model. In this case, the statistical criterion suggests that the model 
should be rejected, and we construct a new and more perfect model. It is often the case 
in practice that Xred > 1 an d the corresponding critical significance level is a < 10%. In 
this case, if we accept the model, we must remember that the obtained parameters and 
their confidence intervals are subject to systematic errors due to either the fact that the 
perfectly exact light curve does not optimally pass through the observed light curve, or the 
presence of a strong correlation in the observational data (note that, in this case, we cannot 
unambiguously judge the systematic error to be due to a correlation in the observational 
data, in contrast to the case when a > 50%). Of course, the obtained model parameters 
and their uncertainties will be less reliable in these cases. The questions raised above are 
considered in more detail elsewhere [Tj [TUl [TU [T3] . 

OBSERVATIONAL DATA 

We analyze here a multicolor transit light curve of the binary system HD 189733, 
which contains an exoplanet [5]. The light curves were obtained with the Hubble Space 
Telescope during three observing runs, on May 22 and 26 and July 14, 2006, each lasting 
five orbital turns of the space observatory. 

The observations of HD 189733 were carried out using the Advanced Camera for 
Surveys (ACS) in the HRC mode. In total, ten light curves with eclipses were obtained in 
the ranges AA = 5500-6000A, 6000-6500A, 6500-7000A, 7000-7500A, 7500-8000A, 
8000 - 8500A, 8500 - 9000A, 9000 - 9500A, 9500 - 10000A, 10000 - 10500A. More 
detailed information about the observational data is available in [5]. When analyzing the 
light curves, we adopted the central wavelengths AA = 5750A, 6250A, 6750A, 7250A, 
7750A, 8250A, 8750A, 9250A, 9750A, 10250A. 

Each light curve contains 675 individual brightness estimates. Figure [T] shows the 
observed light curve in the wavelength range A = 5500 — 6000A (central wavelength 
A = 5750 A). 

We can see from FigJTJthat spots on the surface of the star [5] give rise to considerable 
brightness changes at eclipse phases. There is an appreciable shift in brightness between 
the right and left parts of the light curve due to either a spot on the surface of the star or 
systematic errors in the observations. These features required careful analysis and some 
corrections to the light curve. 

We assumed that the observational uncertainties could be described by a normal dis- 
tribution. The dispersions a 2 of the individual points in the light curve were assumed to 
be the same for all points in the light-curve section used for the fitting. The dispersions 
were determined by averaging the squared differences between the observed points outside 
eclipse and themean brightness outside eclipse, using the left and right branches of the 
light curve both together and separately. We ascribed the resulting values of a as the rms 
deviations of the individual brightness measurements when fitting the corresponding part 
of the light curve. 

We excluded part of the first observing run outside the eclipse at phases 9 < 160°, 
where an obvious systematic shift was observed relative to all other brightness measure- 
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Table 1: rms deviations and the brightness outside eclipse for the light curves 



A(A) 




a both 


M both 




a left 


M left 


^riqht 


^riqht 


M r i a ht 


5750 


0.99998229 


0.00007672 


146 


0.99997273 


0.000078 


79 


0.99999736 


0.000065 


67 


6250 


0.99997042 


0.0003168 


146 


0.99994661 


0.00035 


79 


1.0000079 


0.00026 


67 


6750 


0.99994702 


0.0002085 


145 


0.9999172 


0.00022 


78 


0.99999414 


0.00017 


67 


7250 


0.99996076 


0.0001496 


145 


0.99993531 


0.00016 


78 


1.000001 


0.0001 


67 


7750 


0.99997562 


0.0001421 


146 


0.99996303 


0.00014 


79 


0.99999546 


0.00014 


67 


8250 


0.99999351 


0.0001321 


147 


0.99999245 


0.00013 


80 


0.99999517 


0.00013 


67 


8750 


1.0000042 


0.0001558 


147 


1.0000095 


0.00016 


80 


0.9999958 


0.00015 


67 


9250 


0.99999268 


0.0001766 


147 


0.99999 


0.00018 


80 


0.9999969 


0.00017 


67 


9750 


1.0000026 


0.0002123 


145 


1.0000063 


0.00022 


78 


0.99999678 


0.0002 


67 


10250 


1.0000438 


0.0003353 


148 


1.0000711 


0.00036 


80 


1.0000006 


0.00028 


68 



ments outside the eclipse. 

The mean brightness L outside the eclipse, rms deviations of the individual brightness 
measurements a, and the number of observational points in the eclipse M for each light 
curve (labeled "both" for the left and right branches together, "left" for the left branch 
alone, and "right" for the right branch alone) are given in Table . The accuracy of the 
observed light curve is ~ 10~ 4 , which is ~ 1% of the eclipse depth. 

The main feature of the observed eclipse light curves is the distortion at the minimum 
(at phases 180° < 9 < 190° for the points from the first observing run), due to spots on 
the surface of the star [5]. Moreover, even if this effect is excluded from consideration, 
the light curve remains appreciably distorted near the minimum at 9 ~ 180°. Therefore, 
in addition to fitting using all the points in the light curve, we also analyzed the left 
(descending) branch of the light curve without 9 > 180° and the right (ascending) branch 
without points 9 < 180°. 

The resulting light curves that were analyzed are presented in For ease of viewing, 
the zero points of the light curves for different A are shifted relative to each other. The 
theoretical light curves and the corresponding residual curves calculated using the best-fit 
model with quadratic limb-darkening are shown in the same figure. 

FITTING OF THE LIGHT CURVES USING THE LINEAR LIMB- 
DARKENING LAW 

The parameters derived from the observed light curves are the radius of the star 
r s , radius of the exoplanet r p , orbital inclination i, and coefficient x in the linear limb- 
darkening law. The HD 189733 system contains a star of spectral type K2V [7J. The 
orbital period was taken to be P or h = 2 d . 218581 [5], the planet-to-star mass ratio to 
be q = m p /m s = 0.014 j5], the orbit of the system to be circular, and the radius of the 
relative orbit to be unity. The residual was minimized simultaneously over all parameters. 
Observational points from the phase interval covering only the eclipsed part of the light 
curve were used (the numbers of such points M for corresponding parts of the light curves 
are given in Table ^. The parts of the light curves outside the eclipse were not used in 
the fitting, because they were used to independently determine the rms deviations of the 
points in the light curves (Table,). 
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Differential- Correction Method 

The results of fitting of the left and right parts of the transit light curve and the light 
curve as a whole are given in Tables [2} [31 and HI respectively. The central parameter 
values and their 2a uncertainties were obtained using the differential-correction method. 
Since we used a model that most often proved to be "bad" (see below), we must provide 
2a parameter uncertainties, corresponding to 7 = 95.5%. Moreover, when we determine 
the parameter uncertainties using the differential-correctionmethod in a multiparametric 
model, the true value of each parameter falls into the corresponding error interval with 
the given probability regardless of whether the true values of the other parameters fall 
into their error intervals. 

The question thus arises: what is the probability of all the parameters falling into 
their error intervals? It is natural that we want this probability to be at least at the 
specified confidence level. If we consider the combined central parameter values to be 
a P-dimensional random value varying in the P-dimensional parameter space 0\ . . . 0p , 
then, when the true value 0i of one parameter 0i falls into the corresponding error interval, 
the point associated with the combined true values will fall into the area between the 
two (P-l)-dimensional hyperplanes 0i = 0f — nai and 0i = 0f + KO"j (k is a coefficient 
corresponding to the chosen confidence level). 

Therefore, the differential-correctionmethod is essentially equivalent to finding a P- 
dimensional confidence area of this type. When the true parameter values fall into their 
corresponding error intervals, the point associated with the combined true parameter 
values falls into the P-dimensional parallelepiped centered on the point 0\,...,0p and 
with sides of length 2nai, . . . , 2nap. It is obvious that the probability for the combined 
true values to fall into this parallelepiped will be less than the probability for the true 
parameter values to independently fall into their error intervals, that is, less than the given 
probability. To ensure that a point in the P-dimensional parameter space corresponding to 
the combined true parameter values falls into the P-dimensional parallelepiped centered on 
the point 01, . . . , 0p and with sides proportional to a±, . . . ,ap with the given probability 
7 (i.e., to achieve a situation where the true parameter values all fall simultaneously into 
their corresponding intervals with the given probability 7), we must choose the intervals 
(dimensions of the P-dimensional parallelepiped) to be k^a,, where k > 1. 

The probability for a P-dimensional random point to fall into the P-dimensional par- 
allelepiped is equal to the P-fold integral of the multidimensional density distribution 
over the specified parallelepiped. It is not only the rms deviations of the random values 
01, . . . ,0p that we need to know before defining the type of density function (and deter- 
mining the integral of the function). If a model is a linear function of the parameters 
01, ... , j5p (which have the Gaussiandistributed central values 0\, . . . , ftp) [10]), this dis- 
tribution function will be written where K is a normalizing coefficient and the matrix A 
is the inverse of the covariance matrix cov(/3 i , 0-) 

1 Thc exponent in is the difference between the residual functional as a function of the parameters 
and the minimum value of this functional |10j . and each level surface of the function f (which is a P- 
dimensional ellipsoid) bounds a confidence area obtainedwith a certain confidence level using the Xp 
distribution. 
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p 

m,...^ P ) = Ke^= 1 (3) 

Therefore, even if we use a linear model (or a linear approximation) to determine 
the probability for the true parameter values to jointly fall into their error intervals (or 
if we wish to solve for the coefficient k), we must know the covariances of the central 
values. The same values of the rms deviations o~i, . . . , <jp can correspond to different joint 
probabilities (different coefficients k) in different cases. In any case, however, k/tcxj will 
be less than the corresponding projections of the confidence area obtained using the Xp 
distribution for the given tS 

For example, for 7 = 0.955.. (that is, with k = 2), the projection of the confidence 
area obtained using the xp distribution is Ap = 3.12<r for a four-parameter model and 
Ap = 3.36a for a five-parameter model [TO]. Hence, at the "2cr" (7 = 0.955..) confidence 
level, k < 1.56 for a four-parameter linear model and k < 1.68 for a fiveparameter linear 
model. Note that k < 2.17 for a fourparameter model and k < 2.42 for a five-parameter 
model at the 1 confidence level (7 = 0.68). 

In addition to the parameter values, tablefSJ tableJHl tableJUalso contain the minimum 
reduced x 2 values x 2 e d = (Xjw-p)mW (M — P), and the critical significance levels a . As 
was already noted, the Xred values can be used to test whether the model adequately fits 
the observational data. This indicates (Tables [2j [3], H|) that our model was "bad" in most 
cases. When fitting the right branch of the light curve (Table [3]). 

Fitting of the left branch of the light curve (Table [2]) yielded Xred < 1 f° r three of 
ten wavelengths, with the corresponding ao > 50%. This suggests that systematic errors 
strongly affect the left branch of the light curve. As was noted above, it is these difficulties 
({xled < 1 an d xled) with fitting the light curves of the HD 189733 system that force us 
to adopt the confidence level 7 = 95.5% (which corresponds to 2a in the differential- 
correction method) rather than 68%, indicative of "good" models. For the same reason, 
we were not able to construct "exact" confidence areas D using the x\a distribution and 
to provide the most conservative estimates of the "external" parameter uncertainties for 
the right branch of the light curve (Table [3]) and for the entire light curve (Table H]). We 
restricted our consideration to producing asymptotic confidence areas for the parameters 
using the Xp statistic. We were able to construct the confidence areas for the left branch 
of the light curve using both the Xp an d the x\i distributions. 

Confidence- Area Method 

We used the xp an d x\i distributions for the confidence-area method. We took 7 = 
95.5%, which corresponds to 2<r for the differential- correct ion method, where o is the rms 
deviation. The Xm statistic was used only for the left branch of the light curve. The 
fitting results for the left branch of the light curve obtained using the confidence-area 
method based on the Xm statistic are given in Table |HJ 

1 Since these projections are dimensions of the P-dimensional parallelepiped that are proportional to 
at, ... ^ ap, the probability for them to encompass the true parameter values is greater than 7. 
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Table 2: Fitting of the left branches of the observed light curves for HD 189733 [6] with 
the linear limb-darkening law. (Parameter uncertainties estimated using the differential- 
correction method are given at the 2a level. The two last columns give the reduced x 2 
and the corresponding oq.) 



x(A) 




2<Test (fj) 


V 


2<r eB t 




2<r OB t (* c )(°) 


x c 


2(7 cs t (X C ) 


x 2 red 


a 


5750 


0.11184 


0.00054 


0.01760 


0.00011 


85.715 


0.054 


0.555 


0.017 


2.3457 


1.9 • 10~ u 


6250 


0.1130 


0.0011 


0.01790 


0.00025 


85.60 


0.11 


0.609 


0.037 


0.51644 


0.99996 


6750 


0.11185 


0.00090 


0.01765 


0.00020 


85.697 


0.091 


0.615 


0.028 


0.78673 


0.956 


7250 


0.11133 


0.00073 


0.01751 


0.00016 


85.754 


0.073 


0.590 


0.022 


0.96149 


0.695 


7750 


0.11172 


0.00066 


0.01758 


0.00014 


85.716 


0.065 


0.556 


0.020 


1.1113 


0.345 


8250 


0.11219 


0.00073 


0.01768 


0.00015 


85.693 


0.072 


0.540 


0.023 


1.5651 


0.0031 


8750 


0.11225 


0.00072 


0.01765 


0.00015 


85.682 


0.071 


0.519 


0.023 


1.0306 


0.53 


9250 


0.11186 


0.00080 


0.01756 


0.00016 


85.731 


0.079 


0.490 


0.025 


1.0506 


0.48 


9750 


0.1110 


0.0010 


0.01740 


0.00021 


85.79 


0.10 


0.479 


0.033 


1.1975 


0.19 


10250 


0.1117 


0.0013 


0.01755 


0.00026 


85.75 


0.13 


0.497 


0.041 


0.69904 


0.99 



Table 3: Fitting of the right branches of the observed light curves for HD 189733 [6] with 
the linear limb-darkening law. (Parameter uncertainties estimated using the differential- 
correction method are given at the 2a level. The two last columns give the reduced \ 2 
and the corresponding a .) 



A(A) 




2<r ost (r c s ) 


V 


2<W 


i c (°) 


2 CTrat (i c )(°) 


x c 


2<T ost (x a ) 


xl ed 


a 


5750 


0.11159 


0.00069 


0.01738 


0.00015 


85.793 


0.072 


0.578 


0.023 


4.2535 





6250 


0.1124 


0.0015 


0.01761 


0.00034 


85.71 


0.16 


0.590 


0.053 


1.3380 


0.075155 


6750 


0.1118 


0.0012 


0.01746 


0.00028 


85.77 


0.13 


0.631 


0.041 


1.9672 


0.000029131 


7250 


0.11164 


0.00085 


0.01741 


0.00019 


85.791 


0.091 


0.617 


0.028 


2.6644 


1.2811 ■ 10~ 10 


7750 


0.11173 


0.00086 


0.01740 


0.00019 


85.784 


0.091 


0.593 


0.029 


1.4333 


0.030473 


8250 


0.11177 


0.00078 


0.01742 


0.00017 


85.786 


0.082 


0.570 


0.026 


1.3065 


0.098526 


8750 


0.11185 


0.00088 


0.01738 


0.00019 


85.771 


0.092 


0.534 


0.030 


1.4180 


0.035518 


9250 


0.11104 


0.00092 


0.01725 


0.00019 


85.845 


0.097 


0.521 


0.032 


1.1435 


0.31474 


9750 


0.1108 


0.0012 


0.01727 


0.00025 


85.84 


0.12 


0.497 


0.042 


1.4684 


0.021212 


10250 


0.1110 


0.0013 


0.01726 


0.00027 


85.83 


0.13 


0.506 


0.045 


0.84083 


0.89432 



Table 4: Joint fitting of the right and left branches of the observed light curves for HD 
189733 [6] with the linear limbdarkening law. (Parameter uncertainties estimated using 
the differential-correctionmethod are given at the 2a level. The two last columns give the 
reduced x 2 an d the corresponding a )- 



A(A) 




2o- eB t (rj) 


V 


2*ast (r|) 


; c o 


2<r c5t 


x c 


2<r<,st (x a ) 


x 2 rcd 


"0 


5750 


0.11192 


0.00059 


0.01755 


0.00012 


85.726 


0.060 


0.556 


0.019 


5.5168 





6250 


0.1129 


0.0011 


0.01780 


0.00024 


85.64 


0.11 


0.600 


0.036 


1.0071 


0.55 


6750 


0.11204 


0.00086 


0.01761 


0.00019 


85.706 


0.088 


0.614 


0.028 


1.5321 


0.00010 


7250 


0.11171 


0.00068 


0.01753 


0.00015 


85.740 


0.070 


0.591 


0.022 


1.9306 


6 ■ 10- 10 


7750 


0.11196 


0.00068 


0.01756 


0.00014 


85.717 


0.069 


0.561 


0.022 


2.0916 


lO" 12 


8250 


0.11215 


0.00066 


0.01760 


0.00014 


85.715 


0.067 


0.545 


0.022 


2.2529 


4 ■ 10~ 15 


8750 


0.11230 


0.00073 


0.01758 


0.00015 


85.695 


0.073 


0.515 


0.024 


2.0264 


2 ■ 10 -11 


9250 


0.11166 


0.00070 


0.01747 


0.00014 


85.758 


0.071 


0.494 


0.023 


1.4922 


0.00029 


9750 


0.11103 


0.00080 


0.01736 


0.00016 


85.804 


0.080 


0.482 


0.026 


1.3541 


0.0063 


10250 


0.11150 


0.00096 


0.01745 


0.00020 


85.765 


0.097 


0.493 


0.032 


0.79637 


0.98 



The fitting results obtained using the confidencearea method based on the Xp distri- 
bution are given in TablesEl El [71 Here, the best-fit values of r p , r s , i, x are given together 
with their uncertainties, characterized by the projections of the asymptotic confidence area 
D in the four-parameter space onto the (the confidence intervals). To 

facilitate comparison with the uncertainties obtained in the differential-correction method 
(Tables 0- HI), Tables O - [3 show the values Ap, equal to half the confidence intervals. 
The tables for the confidencearea method present parameter values corresponding to the 
middles of the confidence intervals, rather than with the central parameter values. The 
probability that the exact value of each parameter r p , r s , i, x is encompassed within the 
corresponding confidence interval exceeds 95.5%. The probability that the exact solution 
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(the combination of all parameters r p , r s , i, x) is encompassed within the asymptotic con- 
fidence area is close to the given probability 7 = 95.5%, since the number of points in the 
light curve is large (M > 70). The probability that the exact solution is simultaneously 
encompassed within all confidence intervals exceeds the given 7 = 95.5%. 

Table 5: Fitting of the left branches of the observed light curves for HD 189733 [6] 
with the linear limb- darkening law. (Parameter uncertainties were estimated using the 
confidence-area method based on the Xp distribution, with 7 = 0.955.) 



A(A) 


r s 


A P (r s ) 


r p 


Ap (r p ) 


i 


Ap (i) 


X 


Ap (x) 


5750 


0.11184 


0.00055 


0.01760 


0.00011 


85.715° 


0.054° 


0.555 


0.017 


6250 


0.1130 


0.0025 


0.01788 


0.00055 


85.61° 


0.25° 


0.607 


0.080 


6750 


0.1118 


0.0016 


0.01764 


0.00035 


85.70° 


0.16° 


0.614 


0.049 


7250 


0.1113 


0.0012 


0.01751 


0.00025 


85.76° 


0.12° 


0.590 


0.035 


7750 


0.11171 


0.00094 


0.01758 


0.00020 


85.717° 


0.092° 


0.556 


0.029 


8250 


0.11217 


0.00091 


0.01767 


0.00019 


85.695° 


0.090° 


0.540 


0.029 


8750 


0.1123 


0.0011 


0.01765 


0.00022 


85.68° 


0.10° 


0.518 


0.034 


9250 


0.1119 


0.0012 


0.01756 


0.00025 


85.73° 


0.12° 


0.490 


0.040 


9750 


0.1110 


0.0015 


0.01739 


0.00029 


85.79° 


0.14° 


0.477 


0.046 


10250 


0.1117 


0.0026 


0.01754 


0.00053 


85.75° 


0.26° 


0.493 


0.082 



Table 6: Fitting of the right branches of the observed light curves for HD 189733 [6] 
with the linear limb- darkening law. (Parameter uncertainties were estimated using the 
confidence-area method based on the Xp distribution, with 7 = 0.955.) 



A(l) 


r s 


Ap (r s ) 


r P 


Ap (r p ) 


i 


Ap (i) 


X 


Ap (x) 


5750 


0.11159 


0.00052 


0.01738 


0.00011 


85.794° 


0.055° 


0.578 


0.018 


6250 


0.1124 


0.0019 


0.01760 


0.00043 


85.72° 


0.20° 


0.589 


0.069 


6750 


0.1118 


0.0013 


0.01745 


0.00031 


85.77° 


0.14° 


0.631 


0.046 


7250 


0.11164 


0.00082 


0.01741 


0.00019 


85.792° 


0.089° 


0.617 


0.028 


7750 


0.1117 


0.0011 


0.01739 


0.00025 


85.79° 


0.12° 


0.593 


0.038 


8250 


0.1118 


0.0010 


0.01742 


0.00023 


85.79° 


0.11° 


0.570 


0.036 


8750 


0.1118 


0.0011 


0.01737 


0.00024 


85.77° 


0.12° 


0.533 


0.039 


9250 


0.1110 


0.0013 


0.01725 


0.00028 


85.85° 


0.14° 


0.520 


0.046 


9750 


0.1108 


0.0015 


0.01726 


0.00031 


85.85° 


0.16° 


0.496 


0.053 


10250 


0.1110 


0.0022 


0.01725 


0.00045 


85.83° 


0.23° 


0.503 


0.076 



THE PLANET'S RADIUS AND LINEAR LIMB-DARKENING COEFFI- 
CIENT AS FUNCTIONS OF WAVELENGTH IN THE LINEAR LAW 

Our model for the system HD 189733 with a linear limb-darkening law proved to be 
"bad" for most wavelengths. This is not surprising, since the light curves trace spots 
on the surface of the star j6], which are not taken into account in the model. Only for 
the left branch of the light curve can ourmodel with linear limb- darkening be rejected 
for most wavelengths at a fairly high significance level, so that we can estimate the most 
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Table 7: Joint fitting of the right and left branches of the observed light curves for HD 
189733 [6] with the linear limb darkening law. (Parameter uncertainties were estimated 
using the confidence-area method based on the x 2 p distribution, with 7 = 0.955). 



x(A) 


1 s 


Ap fr e ) 


V 


Ap (r n ) 


i 


Ap (i) 


X 


Ap (x) 

— 1 \ J 


5750 


0.11192 


00039 


017552 


000083 


85.726° 


0.040° 


0.556 


0.013 


6250 


0.1129 


0.0017 


0.01780 


0.00036 


85.64° 


0.17° 


0.598 


0.055 


6750 


0.1120 


0.0011 


0.01761 


0.00024 


85.71° 


0.11° 


0.614 


0.035 


7250 


0.11170 


0.00078 


0.01752 


0.00017 


85.743° 


0.080° 


0.591 


0.025 


7750 


0.11194 


0.00072 


0.01755 


0.00015 


85.719° 


0.072° 


0.562 


0.023 


8250 


0.11215 


0.00068 


0.01760 


0.00014 


85.715° 


0.069° 


0.545 


0.022 


8750 


0.11229 


0.00078 


0.01758 


0.00016 


85.698° 


0.078° 


0.515 


0.026 


9250 


0.11166 


0.00093 


0.01747 


0.00019 


85.759° 


0.095° 


0.493 


0.031 


9750 


0.1110 


0.0010 


0.01736 


0.00021 


85.81° 


0.10° 


0.482 


0.034 


10250 


0.1115 


0.0017 


0.01745 


0.00036 


85.77° 


0.18° 


0.491 


0.058 



Table 8: Fitting of the left branches of the observed light curves for HD 189733 [6J with the 
linear limb-darkening law. (Parameter uncertainties were estimated using the confidence- 
area method based on the Xm distribution, with 7 = 0.955). 



A(A) 


r s 


A M (r s ) 


r p 


Am (r p ) 


i 


A M (i) 


X 


Am (x) 


6250 


0.1127 


0.0066 


0.0178 


0.0014 


85.67° 


0.67° 


0.59 


0.22 


6750 


0.1118 


0.0034 


0.01761 


0.00074 


85.71° 


0.34° 


0.61 


0.10 


7250 


0.1113 


0.0020 


0.01750 


0.00043 


85.76° 


0.20° 


0.589 


0.061 


7750 


0.1117 


0.0013 


0.01758 


0.00027 


85.72° 


0.13° 


0.556 


0.040 


8750 


0.1122 


0.0017 


0.01765 


0.00035 


85.69° 


0.17° 


0.517 


0.055 


9250 


0.1118 


0.0019 


0.01755 


0.00039 


85.74° 


0.19° 


0.488 


0.061 


9750 


0.1110 


0.0016 


0.01739 


0.00032 


85.79° 


0.16° 


0.477 


0.051 


10250 


0.1116 


0.0059 


0.0175 


0.0012 


85.78° 


0.61° 


0.47 


0.19 



conservative "external" parameter uncertainties for 7 = 95.5% using the Xm distribution 
(TableEJ). 

The exoplanet-to-star radius ratio r p /r s in the HD 189733 system was obtained in [6] 
as a function of wavelength. Our results based on our analysis of the complete light curves 
are given in Tabled together with the results obtained in [6|. Since we are searching for 
the individual parameters r p and r s , rather than their ratio r p /r s , the uncertainty in 
r p /r s was obtained using the differential-correctionmethod as an uncertainty for a new 
parameter. The general formula for estimating the rms deviations when changing variables 
is given in [I]. In our case, 

»wu = y / (^) 2 ^>- 2 «-('vs> (JJ) (|) + 

where cov(-, •) is the procedure for searching for the covariances of random values. 
Our uncertainties are larger than those of j6] , where fits were obtained for the geometrical 
parameters r s ,r p ,i, with the limb- darkening coefficients being fixed. Since we fit the 
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light curves allowing both the geometrical parameters and the limb-darkening coefficients 
to vary, the larger number of degrees of freedom in our model yields larger parameter 
uncertainties. 

Figure [3] indicates a certain increase in the planet's radius with decreasing wavelength. 
This effect is significant at the la level, but not the 2a level. Our values of the exoplanet 
radius are systematically larger than those determined in [5] (by 0.3%); this is due to the 
different normalization procedures applied in the analysis of the light curve and the fact 
that we searched for the limb-darkening coefficient together with other parameters, rather 
than fixing this coefficient. 

Table 9: Comparison of the r p /r s values as a function of wavelength obtained in this 
paper and in [BJ. (The la uncertainties are given). 



A(A) 




m 


5750 


0.156894 ± 0.000281 


0.156903 ±0.000095 


6250 


0.157762 ± 0.000658 


0.156744 ±0.000065 


6750 


0.157387 ± 0.000458 


0.156552 ±0.000057 


7250 


0.157045 ± 0.000341 


0.156388 ±0.000059 


7750 


0.156917 ± 0.000358 


0.156501 ±0.000064 


8250 


0.156984 ± 0.000348 


0.156210 ±0.000073 


8750 


0.156569 ± 0.000388 


0.156147 ±0.000081 


9250 


0.156451 ± 0.000391 


0.156120 ±0.000092 


9750 


0.156361 ± 0.000463 


0.156097 ±0.000125 


10250 


0.156389 ± 0.000611 


0.155716 ±0.000218 



Table 10: Fitting of the left branches of the observed light curves for HD 189733 [6J 
with the quadratic limb-darkening law. (Parameter uncertainties estimated using the 
differential- correct ion method are given at the 2a level. The two last columns present the 
values of the reduced \ 2 an d the corresponding critical significance levels ao). 



x(A) 




2o- oat 


r p 


2<7™t (rg) 




2<r OB t (i c )(°) 


xl 


2tr es t («i) 


vl 


2<r OBt (yf) 


x 2 red 


a 


5750 


0.11207 


0.00083 


0.01754 


0.00022 


85.721 


0.077 


0.44 


0.16 


0.18 


0.26 


2.3115 


8 ■ 10~ a 


6250 


0.1132 


0.0020 


0.01782 


0.00050 


85.62 


0.17 


0.50 


0.37 


0.16 


0.60 


0.52219 


0.99997 


6750 


0.1120 


0.0014 


0.01759 


0.00035 


85.70 


0.12 


0.52 


0.26 


0.14 


0.42 


0.79144 


0.96 


7250 


0.11144 


0.00094 


0.01748 


0.00024 


85.758 


0.088 


0.53 


0.20 


0.08 


0.32 


0.97038 


0.70 


7750 


0.11189 


0.00091 


0.01754 


0.00023 


85.720 


0.083 


0.47 


0.20 


0.13 


0.30 


1.1150 


0.37 


8250 


0.1125 


0.0022 


0.01754 


0.00048 


85.72 


0.16 


0.32 


0.24 


0.33 


0.48 


1.5210 


0.0074 


8750 


0.1123 


0.0010 


0.01757 


0.00026 


85.704 


0.094 


0.42 


0.23 


0.14 


0.35 


1.0401 


0.54 


9250 


0.1121 


0.0013 


0.01748 


0.00036 


85.74 


0.12 


0.34 


0.26 


0.23 


0.41 


1.0466 


0.53 


9750 


0.1112 


0.0013 


0.01737 


0.00034 


85.79 


0.12 


0.41 


0.33 


0.10 


0.49 


1.2103 


0.198 


10250 


0.1131 


0.0044 


0.0175 


0.0012 


85.71 


0.29 


0.05 


0.43 


0.65 


0.45 


0.66148 


0.997 



Note that the difference between the exoplanet radius in HD 209458 obtained in the 
blue and the red does not seem to be significant [lj: — 0.0139 ± 0.0003 for A = 3750 
(the mean value for A = 3201, 3750, 4300AA) and = 0.0138 ± 0.0002 for A = 8732 (the 
mean value for A = 7755, 8732, 9708 AA); the 2a uncertainties are indicated. 

Let us now consider the limb- darkening coefficients in the linear limb-darkening law as 
a function of wavelength. Figures SI El El El show the observed limbdarkening coefficients 
x versus the wavelength A, together with the theoretical linear dependences of the limb- 
darkening coefficients on in the photometric systems ugriz and UBVRIJ j2[ IT5] . The 
observed relation x(X) agrees qualitatively with the theoretical one: the observed x(X) 
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Table 11: Fitting of the right branches of the observed light curves for HD 189733 [6] 
with the quadratic limb-darkening law. (Parameter uncertainties estimated using the 
differential- correct ion method are given at the 2a level. The two last columns present the 
values of the reduced x 2 an d the corresponding critical significance levels oq). 



A(A) 




2<T„ st (rj) 


V 


2<r eB t (r|) 


i c (°) 


2<7 OB t (i c )(°) 




2cr CB t (sf) 


vl 


2CT ost (vi) 




a 


5750 


0.11094 


0.00053 


0.017531 


0.000095 


85.781 


0.044 


0.86 


0.20 


-0.42 


0.27 


3.7576 





6250 


0.1113 


0.0011 


0.01794 


0.00016 


85.681 


0.078 


1.13 


0.44 


-0.78 


0.58 


1.1957 


0.26 


6750 


0.11080 


0.00089 


0.01765 


0.00014 


85.761 


0.070 


1.01 


0.33 


-0.56 


0.45 


1.8492 


2.6 • 10 -4 


7250 


0.11115 


0.00072 


0.01754 


0.00014 


85.772 


0.060 


0.84 


0.25 


-0.33 


0.35 


2.5861 


1.27 ■ 10 -9 


7750 


0.11109 


0.00069 


0.01752 


0.00013 


85.776 


0.058 


0.85 


0.25 


-0.38 


0.35 


1.3392 


0.089 


8250 


0.11076 


0.00050 


0.017623 


0.000080 


85.778 


0.040 


0.99 


0.20 


-0.61 


0.27 


0.99821 


0.65 


8750 


0.11116 


0.00069 


0.01753 


0.00012 


85.761 


0.056 


0.85 


0.27 


-0.46 


0.37 


1.3138 


0.11 


9250 


0.11083 


0.00089 


0.01731 


0.00018 


85.839 


0.079 


0.63 


0.28 


-0.16 


0.40 


1.1495 


0.338 


9750 


0.1120 


0.0024 


0.01735 


0.00061 


85.78 


0.21 


0.31 


0.44 


0.26 


0.73 


1.5716 


0.0089 


10250 


0.1124 


0.0033 


0.01736 


0.00076 


85.74 


0.27 


0.26 


0.48 


0.33 


0.87 


0.88100 


0.86 



values decrease, on average, with increasing wavelength. However, the observed x(X) are 
systematically lower than the theoretical values. This is in qualitative agreement with 
the results obtained for the HD 209458 system [TJ H] . 

However, the behavior of the observed x(X) differs from the x(X) relation obtained 
for HD 209458, where the difference between the observed and theoretical x(X) values 
increases monotonically with A; for the HD 189733 system, the difference is largest at the 
shortest wavelengths and decreases toward longer wavelengts. Figures H] and indicate 
that the results obtained when fitting the left and right branches of the HD 189733 light 
curve separately agree with each other, despite the systematic shift in brightness between 
these light curves. Therefore, we conclude that our results for the limb-darkening of the 
star are stable against systematic errors affecting the light curve in the HD 189733 system. 

Table 12: Joint fitting of the right and left branches of the observed light curves for HD 
189733 [6] with the quadratic limb-darkening law. (Parameter uncertainties estimated 
using the differential- correct ion method are given at the 2a level. The two last columns 
present the values of the reduced x 2 an d the corresponding critical significance levels a ). 



A(A) 




2o- st (^a) 


r p 


2a OBt (rfj 


i c (°) 


2<Tost (i c )(°) 


xl 


2o- OB t (zf) 


vi 


2o- ost (vi) 






5750 


0.11210 


0.00067 


0.01760 


0.00015 


85.705 


0.058 


0.57 


0.18 


-0.02 


0.27 


5.5682 





6250 


0.11228 


0.00094 


0.01792 


0.00017 


85.637 


0.072 


0.85 


0.33 


-0.36 


0.45 


0.99780 


0.608 


6750 


0.11183 


0.00086 


0.01768 


0.00018 


85.697 


0.072 


0.72 


0.26 


-0.15 


0.37 


1.5358 


0.00013 


7250 


0.11156 


0.00071 


0.01756 


0.00015 


85.737 


0.061 


0.66 


0.20 


-0.10 


0.29 


1.9376 


8 ■ lO" 10 


7750 


0.11180 


0.00072 


0.01759 


0.00015 


85.716 


0.061 


0.63 


0.21 


-0.09 


0.30 


2.0992 


2.5 ■ 10~ 12 


8250 


0.11201 


0.00069 


0.01764 


0.00015 


85.712 


0.058 


0.61 


0.20 


-0.10 


0.29 


2.2615 


5 ■ 10~ 15 


8750 


0.11216 


0.00077 


0.01761 


0.00016 


85.695 


0.065 


0.58 


0.24 


-0.10 


0.34 


2.0344 


2.7 ■ 10 -11 


9250 


0.11174 


0.00090 


0.01743 


0.00022 


85.766 


0.084 


0.43 


0.23 


0.09 


0.34 


1.4995 


3 ■ 10" 4 


9750 


0.11105 


0.00091 


0.01736 


0.00021 


85.804 


0.083 


0.47 


0.25 


0.01 


0.37 


1.3637 


0.0063 


10250 


0.1129 


0.0023 


0.01733 


0.00073 


85.74 


0.18 


0.00 


0.35 


0.72 


0.38 


0.77561 


0.98971 



HD 189733 and HD 209458 have qualitatively different functions x(A), probably due 
to the fact that there were spots on the surface of the star in the HD 189733 system at 
the observing epoch [5j |6]. Therefore, additional observations of HD 189733 in periods 
of lower activity would be of great interest. Note that the observed x(X) values in both 
HD 189733 and HD 209458 lie below the theoretical relations, although the observed x(A) 
relations are qualitatively different for them. This result is important for checking modern 
models of thin stellar atmospheres. 
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FITTING OF THE LIGHT CURVES WITH QUADRATIC LIMB DARK- 
ENING 



The results of fitting the light curves using the quadratic limb-darkening law are given 
in Tables [T0l-[T6l The data obtained using the differential-correction method are given in 
Tables [l~0l - fl~2l Here, the central parameter values r p , r s , i, xi, y\ together with their 2a 
uncertainties (7 = 95.5%) are shown. Tables fTUl - [121 also contain the values of xl e d an d a o- 
Fitting of the left and right branches of the light curves (Tables fTUl and ITTT) yields Xred < 1 
and «o > 0.5 for some wavelengths. This indicates that the observed brightnesses may 
be correlated and systematic errors may be present in the observational data. When we 
fit the light curve as a whole (Table 1X21 . our model is "bad" for most wavelengths, and is 
rejected at a very low significance level. The light-curve fits at A = 6250A and A = 10250A 
yield Xred — 1 anc ^ a o > 0.5, suggesting a strong correlation of the observational points. 
All these factors lead us to adopt 2a est rather than a est to estimate the uncertainties of 
the model parameters for the quadratic limb-darkening law, as was done for the linear 
limb-darkening model. 



Table 13: Fitting of the left branches of the observed light curves for HD 189733 [6] 
with the quadratic limb-darkening law. (Parameter uncertainties were estimated using 
the confidence-area method based on the \ 2 p distribution, with 7 = 0.955.) 



A(A) 


r s 


A P (r s ) 


r p 


Ap (r p ) 


i 


Ap (i) 


X\ 


Ap (xi) 


2/1 


Ap (yi) 


5750 


0.11212 


0.00067 


0.01754 


0.00016 


85.722° 


0.061° 


0.43 


0.18 


0.17 


0.27 


6750 


0.1123 


0.0020 


0.01757 


0.00047 


85.71° 


0.18° 


0.51 


0.48 


0.13 


0.73 


7250 


0.1116 


0.0013 


0.01747 


0.00031 


85.76° 


0.12° 


0.51 


0.32 


0.07 


0.49 


7750 


0.1120 


0.0013 


0.01752 


0.00028 


85.72° 


0.11° 


0.45 


0.32 


0.08 


0.48 


8250 


0.1126 


0.0012 


0.01754 


0.00029 


85.72° 


0.12° 


0.32 


0.34 


0.29 


0.51 


8750 


0.1125 


0.0014 


0.01758 


0.00030 


85.69° 


0.12° 


0.41 


0.37 


0.09 


0.54 


9250 


0.1122 


0.0015 


0.01745 


0.00038 


85.76° 


0.15° 


0.32 


0.45 


0.20 


0.67 


9750 


0.1112 


0.0018 


0.01736 


0.00039 


85.80° 


0.16° 


0.43 


0.39 


0.13 


0.62 



Table 14: Fitting of the left branches of the observed light curves for HD 189733 [6] 
with the quadratic limb- darkening law. (Parameter uncertainties were estimated using 
the confidence-area method based on the Xm distribution, with 7 = 0.955). 



X(A) 


r s 


Am (r s ) 


r P 


Am (r p ) 


i 


Am (i) 




Am (xi) 


2/1 


Am (2/1) 


6750 


0.1150 


0.0066 


0.01785 


0.00051 


85.77° 


0.41° 


1.27 


0.81 


1.3 


1.3 


7250 


0.1115 


0.0022 


0.01758 


0.00039 


85.77° 


0.21° 


0.63 


0.40 


0.23 


0.62 


7750 


0.1118 


0.0014 


0.01756 


0.00030 


85.71° 


0.12° 


0.47 


0.27 


0.13 


0.42 


8750 


0.1123 


0.0018 


0.01766 


0.00036 


85.69° 


0.17° 


0.51 


0.36 


0.28 


0.55 


9250 


0.1121 


0.0022 


0.01762 


0.00036 


85.76° 


0.22° 


0.42 


0.35 


0.36 


0.56 


9750 


0.1109 


0.0018 


0.01741 


0.00040 


85.77° 


0.14° 


0.39 


0.40 


0.073 


0.64 



Tables UM HI US] show the results of fitting the light curves of HD 189733 using the 
model with quadratic limb darkening; the parameter uncertainties are estimated using 
the confidence-area method based on the x 2 p distribution. Table [T41 presents the results of 
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Table 15: Fitting of the right branches of the observed light curves for HD 189733 [6] 
with the quadratic limb-darkening law. (Parameter uncertainties were estimated using 
the confidence-area method based on the x 2 p distribution, with 7 = 0.955.) 



X(A) 


' s 


L - y F V s) 


' V 


Ad (r m ) 




Ap (i) 




Ap (xi ) 




Ap (vi ) 


5750 


0.11091 


0.00070 


0.01752 


0.00016 


85.782° 


0.060° 


0.87 


0.18 


-0.44 


0.26 


6250 


0.1109 


0.0028 


0.01784 


0.00064 


85.73° 


0.22° 


1.27 


0.69 


-0.73 


0.97 


6750 


0.1108 


0.0019 


0.01758 


0.00043 


85.77° 


0.16° 


1.02 


0.44 


-0.60 


0.64 


7250 


0.1111 


0.0014 


0.01750 


0.00037 


85.78° 


0.13° 


0.85 


0.37 


-0.39 


0.54 


7750 


0.1111 


0.0013 


0.01751 


0.00030 


85.78° 


0.11° 


0.87 


0.33 


-0.38 


0.48 


8250 


0.1108 


0.0012 


0.01763 


0.00027 


85.778° 


0.098° 


1.01 


0.29 


-0.59 


0.42 


8750 


0.1113 


0.0015 


0.01753 


0.00032 


85.76° 


0.12° 


0.85 


0.38 


-0.45 


0.55 


9250 


0.1109 


0.0016 


0.01731 


0.00041 


85.84° 


0.16° 


0.66 


0.45 


-0.12 


0.67 


9750 


0.1108 


0.0020 


0.01726 


0.00043 


85.85° 


0.18° 


0.55 


0.53 


-0.17 


0.78 



Table 16: Joint fitting of the right and left branches of the observed light curves for 
HD 189733 j6] with the quadratic limbdarkening law. (Parameter uncertainties were esti- 
mated using the confidence-area method based on the Xp distribution, with 7 = 0.955.) 



A(A) 


r s 


A P (r s ) 


r p 


Ap (r„) 


i 


Ap (i) 


X\ 


Ap (xi) 


yi 


Ap (yi) 


5750 


0.11182 


0.00048 


0.01757 


0.00011 


85.724° 


0.042° 


0.60 


0.12 


-0.09 


0.18 


6250 


0.1124 


0.0021 


0.01790 


0.00047 


85.64° 


0.17° 


0.88 


0.53 


-0.35 


0.76 


6750 


0.1119 


0.0014 


0.01765 


0.00031 


85.71° 


0.12° 


0.71 


0.34 


-0.16 


0.50 


7250 


0.11160 


0.00095 


0.01755 


0.00023 


85.741° 


0.085° 


0.65 


0.24 


-0.12 


0.36 


7750 


0.11183 


0.00091 


0.01758 


0.00019 


85.718° 


0.075° 


0.62 


0.23 


-0.10 


0.34 


8250 


0.11204 


0.00085 


0.01763 


0.00019 


85.714° 


0.072° 


0.61 


0.22 


-0.11 


0.32 


8750 


0.11219 


0.00100 


0.01760 


0.00021 


85.697° 


0.082° 


0.58 


0.27 


-0.10 


0.38 


9250 


0.1118 


0.0011 


0.01743 


0.00028 


85.77° 


0.11 


0.45 


0.33 


0.12 


0.48 


9750 


0.1111 


0.0013 


0.01735 


0.00027 


85.81° 


0.11 


0.48 


0.35 


0.01 


0.51 



fitting the left branches of the light curves using the quadratic limb-darkening law. Here, 
the parameter uncertainties were obtained using the confidence-area method based on the 
Xm distribution (7 = 0.955). 

27 AND Vl AS FUNCTIONS OF WAVELENGTH IN THE QUADRATIC 
LIMB-DARKENING LAW 

Figures 181 - [TBI show the observed 37(A) and yi(X) values as functions of wavelength, 
obtained using the differential-correction method. The 2a uncertainties are given (7 = 
95.5%). The theoretical relations 27(A) and yi(X) obtained in [21 US] are also shown 
here. The observed relation 07(A) obtained using the left branch of the light curve agrees 
with the theoretical function 07(A), while the observed xi(A) obtained using the right 
branch of the light curve lies considerably above the theoretical relation (FigsJH] and |9]). 
The observed 27(A) obtained using the whole light curve agrees satisfactorily with the 
theoretical relation within the uncertainties (at the 2a level). 

The observed 2/1 (A) obtained using the left branch of the light curve agrees with the 
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theoretical function, while the observed yi(X) obtained for the right branch of the light 
curve lies considerably below the theoretical relation. The observed yi(X) obtained us- 
ing the whole light curve agrees satisfactorily with the theoretical relation within the 
uncertainties (at the 2a level) (Figs. [TT1- [T3|) . 

Note that the observed relations Xi(X) and yi(X) obtained using the whole light curve 
do not agree with the theoretical relations at the la level: the observed Xi(A) is sys- 
tematically higher than the theoretical relation, and 2/1 (A) is lower than the theoretical 
relation (Figs. ITUl and [TBI . Since our model is formally bad and, moreover, has five rather 
than one parameters, we are forced to take the uncertainties at the 2a level. With these 
uncertainties for the coefficients 21(A) and 2/1 (A), the differences between the observed 
and theoretical 21(A) and 2/1 (A) relations do not seem to be significant. Figures and 
[T5l show the Xi(X) and 2/1 (A) relations obtained when fitting the left and right branches 
of the light curve, with the uncertainties derived using the confidence-area method based 
on the Xp distribution. 

Figures [161 and [T71 illustrate the projections of the confidence areas in the (r s ,r p ) and 
(xi,yi) planes obtained when using the Xp an d Xm distributions. The light curve at 
A9500 ~r~ 10000 A (left branch) was used with the model based on the quadratic limbdark- 
ening law. 

CONCLUSION 

We have analyzed high-accuracy multicolor light curves given in [5], [6] for the star- 
explanet system HD 189733 and determined the radii of the star and planet, the orbital 
inclination, and the coefficients in the linear and quadratic limb-darkening laws across the 
disk of the K2V star. The results of our fitting agree with the data obtained in [51 [6]. 

We have analyzed in detail the limb darkening across the disk of the K2V star based on 
the resulting uncertainties in the coefficients in the linear and quadratic limb-darkening 
laws. We allowed for the presence of spots on the surface of the K2V star [H] by analyzing 
the left and right branches of the light curves separately, as well as the light curve as a 
whole. Moreover, we paid special attention to checking the adequacy of our model and 
verifying how well it agrees with the observational data. Our model proved to be formally 
"bad". When fitting the left and right branches of the light curve separately, our model 
is rejected at a high significance level of a > 50% at some wavelengths; this most likely 
indicates a correlation of the individual observational points in the light curve. When 
fitting the light curve as a whole, our model is rejected at a very low significance level for 
most wavelengths. 

Since our model turned out to be "bad" when applied to HD 189733, we were forced 
to adopt the high confidence level 7 = 95.5% (rather than the 68% relevant for "good" 
models) when estimating the parameter uncertainties. Since our model includes four or 
five parameters instead of only one, we choose the error intervals in the confidence-area 
method such that the probability for them to encompass the true parameter values is 
undoubtedly higher than the given probability (in this case, the given confidence level 
7 = 95.5% is associated with the entire confidence area D rather than only one confi- 
dence interval). Our analysis of the relations 2 (A), £i(A), 2/1 (A) for the limb-darkening 
coefficients deduced from fitting the observations yielded the following results 7 = 95.5%. 
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The observed values of the coefficient x(X) in the linear limb-darkening law for HD 
189733 are systematically below the theoretical relation, with the differences increasing 
with decreasing wavelength A (in contrast to HD 209458, where they increase with in- 
creasing A). The observational coefficients in the quadratic limb-darkening law £i(A) and 
2/i (A) agree satisfactorily within the 2a uncertainties (7 = 95.5%) with the theoretical 
relations developed for one- dimensional thin stellar atmospheres ([2J [3j 115]). 

We emphasize that these conclusions concern the light curves of the HD 189733 system 
when there were spots on the disk of the K2V star. To further invetigate limb darkening 
on this star, additional observational data should be obtained for eclipse curves at epochs 
when the contribution of spots is negligible. 

We confirm the earlier conclusion [5] (at the la level) that the exoplanet radius in- 
creases with decreasing wavelength. This may imply the presence of the atmosphere 
around the exoplanet. 

The authors are grateful to Frederic Pont for providing us with the observational data 
on HD 189733. 
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Fig 1: Light curve 1 of the star HD 189733 in the filter A = 5500 - 6000 taken from ©. 
The light curve has truncated wings and a brightness jump at the minimum (due to a 
spot on the surface of the star). 
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Fig 2: Light curves of HD 189733 at A = 5500 — 10500. The wavelenght increases upward. 
The residuals for the quadratic limb-darkening law with the best-fit parameter values are 
given in the lower part. 
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Fig 3: The planet-to-star radius ratio as a function of wavelength obtained in this paper 
(filled circles) and according to [6] (filled squares). The indicated uncertainties are la. 
Since we do not specify the limb- darkening coefficient, and are searching for it together 
with the other parameters of the problem, our uncertainties are larger. 
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Fig 4: Limb-darkening coefficient x for HD 189733 as a function of wavelength A de- 
rived for the linear limb- darkening law and the left branches of the light curves from [6]. 
The uncertainties of the limb- darkening coefficient were obtained using the differential- 
correction method and are given at the 2a level. The theoretical values of the limb- 
darkening coefficients in the ugriz and UBVRIJ photometric systems were taken from 
Claret2004,Claretl998,Claret2000. 
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Fig 5: Limb- darkening coefficient x for HD 189733 as a function of wavelength A derived 
for the linear limb-darkening law and the right branches of the light curves from [6] . The 
uncertainties of limb-darkening coefficient were obtained using the differential-correction 
method and are given at the 2a level. The theoretical values of the limb-darkening 
coefficient in the ugriz and UB V RIJ photometric systems were taken from [21 [HI [15] . 
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Fig 6: Limb-darkening coefficient x forHD189733 as a function of wavelength A derived 
for the linear limb-darkening law and the joint analysis of the right and left branches of the 
light curves from |6j. The uncertainties of limb-darkening coefficient were obtained using 
the differential-correction method and are given at the 2 level. The theoretical values of 
the limb- darkening coefficient in the ugriz and UBVRIJ photometric systems were taken 
from [2JHUH5]. 
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Fig 7: Same function as in Fig. 6 but with the uncertainties of the limb-darkening co- 
efficients obtained using the confidencearea method based on the x 2 p distribution; the 
uncertainties are given for 7 = 0.955. 
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Fig 8: Same function as in Fig.4 but for the linear coefficient X\ and the quadratic limb- 
darkening law. 
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Fig 9: Same function as in Fig. 8 but obtained from the analysis of the right branches of 
the light curves. 
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Fig 10: Same function as in Fig. 8 but obtained from the joint analysis of the left and 
right branches of the light curves. 
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Fig 11: Same function as in Fig. 8 but for the quadratic coefficient y\. 
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Fig 12: Same function as in Fig. 11 but obtained for the analysis of the right branches of 
the light curves. 
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Fig 13: Same function as in Fig. 11 but obtained for the joint analysis of the left and 
right branches of the light curves. 
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Fig 14: Same function as in Fig. 10 but with the uncertainties in limb-darkening coefficients 
calculated using the confidencearea method based on the Xp distribution; uncertainties 
correspond to 7 = 0.955. 



34 




Fig 15: Same function as in Fig. 14 but for the quadratic coefficient y\. 
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Fig 16: Projections of the confidence area (for 7 = 0.955) obtained using the Xp (dashed 
line) and Xm ( s °hd line) distributions in the plane of the parameters r s ,r p when fitting 
the left branch of the light curve for 9500 — 10000A with the quadratic limb-darkening 
law. The sides of the rectangle correspond to the 2a error intervals obtained using the 
differential-correction method. 
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Fig 17: Projections of the confidence area (for 7 = 0.955) obtained using the Xp (dashed 
line) and xli ( s °hd line) distributions in the plane of the parameters x 1 ,t/i when fitting 
the left branch of the light curve for 9500 — 10000A with the quadratic limb-darkening 
law. The sides of the rectangle correspond to the 2a error intervals obtained using the 
differential-correction method. 
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